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We consider Saffman- Taylor channel flow without surface 
tension on a high-pressure driven interface, but modify the 
usual infinite-fluid in infinite-channel configuration. Here we 
include the treatment of efflux by considering a finite con- 
nected body of fluid in an arbitrarily long channel, with its 
second free interface the efflux of this configuration. We show 
that there is a uniquely determined translating solution for 
the driven interface, which is exactly the 1/2 width S-T so- 
lution, following from correct symmetry for a finite channel 
flow. We establish that there exist no perturbations about 
this solution corresponding to a finger propagating with any 
other width: Selection is unique and isolated. The stability of 
this solution is anomalous, in that all freely impressible per- 
turbations are stabilities, while unstable modes request power 
proportional to their strength from the external agencies that 
drive the flow, and so, in principle, are experimentally control- 
lable. This is very different from the behavior of the usual infi- 
nite fluid: The limit of infinite length is singular, and not that 
of the literature. We argue that surface tension on the efflux 
interface modifies channel- width A according to 1 — 2 A = cr/v 
(i.e. (2ty) 2 B of the literature) with v the velocity of the high- 
pressure tip, but a the surface tension of the efflux. That is, A 
is decreased below 1/2 by the long-range effect of smoothing 
an arbitrarily distant efflux. The perturbation theory created 
here to deal with transport between two free boundaries is 
novel and dependent upon a remarkable symmetry. 



1. INTRODUCTION 

Objects arising from flows, aggregation and so forth 
often have distinctive shapes. (Think, for example, of 
snowflakes.) Upon theoretical study, (the identification 
of all forces and processes at work, and the determination 
of all boundary conditions), the solution will be that of 
this particular shape. In other cases, however, the usual 
form is just one of many possible solutions that can be 
found. The question of what then distinguishes and de- 
termines only the usual, observed form is termed a prob- 
lem of "pattern selection" . Clearly something has been 
left out of the theoretical treatment, which means that 
after a careful treatment, something apt to be subtle re- 
mains to be dealt with. "Left out" itself is subtle, since 
the theoretical considerations, and their solutions, entail 
approximations, and the "truth" can have been left out 
by the approximations. We consider such a problem here. 

Perhaps the earliest studied selection problem arose 
in the 1958 work of Saffman- Taylor (S-T) Here the 
question is the shape of a bubble penetrating into a body 



of viscous fluid. Loosely, a soda-straw is filled with a 
thick fluid which is then forcibly blown out. How has 
the breath of air propagated down the straw, to drive 
out what amount of the fluid? Performed carefully, the 
straw of length Z is flattened down all along its length 
into a long rectangle of width w and length Z, where 2w 
is just about the circumference of the cylindrical straw, 
and a channel of thickness b <C w all that remains of the 
third dimension of the straw's volume. The material of 
the "straw" is transparent glass for visualization. That 
is, there are two w by Z plates of glass, one a gap b above 
the other, and the long edges are sealed. View the w by 
Z rectangles from above with the long length horizontal 
(the x-axis) and the width w along the y-axis. A con- 
trolled source of high pressure air is now applied to the 
left end (say, at x=0), and the driven viscous fluid flows 
out the right end into atmosphere (the efflux). The bub- 
ble of air advancing from the left into the viscous fluid, 
the immiscible high-pressure interface between the driv- 
ing air and viscous fluid, hereafter referred to as H, is 
the object of study. (See Fig.l, but imagine the viscous 
fluid to continue beyond L.) 




FIG. 1. The channel in (x, y) coordinates. The body of 
viscous fluid is shaded. It is bounded between a high-pressure 
interface H with "air" and another, L at low-pressure. 

Invariably, S-T observed the bubble to develop into a 
very long "finger" , that is, a rounded nose furthest to the 
right with long horizontal sides parallel to the long sides 
of the channel trailing back towards the high pressure 
end. This finger is almost perfectly centered within the 
channel, that is, with center-line symmetry. Moreover, 
S-T discovered that with higher and higher pressures, 
and hence increasingly faster bubbles, the width of the 
horizontal sides of the finger approached 1/2 the width 
w of the channel. This is the selected pattern for this 
problem. 

S-T then theoretically determined H of finger shape, 
but of width Aw with A anywhere from to 1, rather 
than just 1 /2. They had no argument as for why only the 
observed value should occur, and so exposed a problem 
of pattern selection. 

Let us be more careful and say more about the exper- 
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iment. There is no reason to use just air to drive the 
viscous fluid. Indeed, any fluid immiscible with the fluid 
initially filling the channel may be employed to drive it. 
The problem is simplest if the driving fluid, like air, has 
a negligible viscosity, so that pressure is spatially uni- 
form throughout it, and hence all along the interface H. 
However, there is always some surface tension along H, 
so that when H is curved, just within the viscous fluid 
along H, pressure is not uniform. It is straightforward to 
see (say on dimensional grounds) that this surface tension 
a can modify the problem only in the combination o~/v 
with v a velocity, say the velocity of the tip of the finger. 
It then follows that the selection of 1/2, occurring for 
high velocities, is selection in the limit a/v — > + . This 
then sets up the question of how this singular limit ac- 
counts for the observed selection. It is singular because, 
as determined in the 1958 paper, with a = 0, there is no 
selection at all. 

There is an extensive literature, culminating in the 
mid-1980s and employing "beyond all orders" of per- 
turbation methodology to verify this singular limit, and 
hence the critical role played by even the slightest of sur- 
face tensions acting along the boundary of the selected 
pattern. It is now the accepted wisdom that all these 
indeterminate selection problems achieve their theoreti- 
cal resolution through additional "constituentive" rela- 
tionships (such as a for S-T) imposed directly upon the 
shape to be so selected. 

Let us again be more careful, and say more about the 
theory. The gap b is viewed as so small, that with a usual 
limit of the viscous fluid equations, this "film" of fluid 
can be treated purely within 2-D. This is a theoretically 
compelling circumstance, since the problem then rapidly 
is noticed to be fully amenable to complex, conformal 
mathematics. [^,^4| Apart from some studies as to how 
finite values of b can modify the 2-D solution, the entire 
literature on S-T flow is 2-D conformal. 

The theoretical and experimental literature neverthe- 
less diverge in an obvious way. While discussing sur- 
face tension on the bubble interface, the experimental 
literature is relatively vague about how the viscous fluid 
emerges from the apparatus at low pressure. It is not 
that this is totally unimportant - care has to be paid to 
terminations, and the reliable reproducibility of fast long 
fingers emerges only after some adjustment of the efflux 
termination. Perhaps it is largely undocumented in con- 
sequence of an implicit belief that far enough away from 
the observed interface, these details should significantly 
have decoupled. 

The theoretical literature has gone much further. The 
entire literature considers only I = oo, namely, an in- 
finitely long channel all filled with the viscous fluid. This 
is evidently a very different geometry from that of the 
experiments that determine selection. Granted that sub- 
tleties can potentially lead to selection or its absence, 
one might wonder that this approximation to very large 
length has an impact upon selection. This question has 
never been asked in the literature. We do so here, and 



indeed discover that the limit I — > oo is also singular. 

Since conformal mathematics is so powerful, we con- 
sider a finite termination of geometrically simplest form. 
Namely, we allow the channel containing the fluid to be 
arbitrarily long, while the viscous fluid, always within it, 
is of a fixed finite volume. Abutting the fluid at its low 
pressure, right end, we again have air; that is, a fluid with 
negligible viscosity, as is the high pressure driving fluid. 
Hence, we now have two interfaces; the high-pressure, 
experimentally observed one H, and the low-pressure, 
efflux one at the right, which we now term L. Both in- 
terfaces are dynamically free surfaces, although they are 
evidently highly coupled, since between them resides the 
connected body of our incompressible viscous fluid. (See 
Fig.l) 

In this finite configuration, the pressure drop across 
the viscous fluid becomes physically finite, whereas it is 
simply infinite in the usual treatment. There is also a 
net flux of fluid. In the infinite treatment it freely can 
be set at any value with impunity, and is conventionally 
set at the constant 1. In our finite case, the impedance 
of the flow, the ratio of net pressure drop by flux is again 
finite and large, as opposed to simply infinite in the usual 
treatment. More importantly, this impedance generally 
changes in time, since the closest distance between the 
two interfaces dynamically changes in time, and gener- 
ally decreases as the interface H grows more bent into 
a finger. In experiment Ref. ||, the time derivative of 
impedance is indeed significantly negative. Although 
impedance variations are buried in an infinite value in 
the usual infinite treatment, it is troublesome that there 
is no trace of its derivative in that literature. 

The reason we have considered a finite body of fluid 
with definite pressures along its free boundaries is that 
there is a mathematical clue to anticipate an influence 
upon shape. The theoretical works on the infinite fluid 
all consider periodic boundary conditions for the cross- 
channel (y) behavior. The channel has two long sides 
a distance w apart, one edge at y=0 and another at 
y=w. Interface H (the only one for the infinite prob- 
lem) has viscous fluid everywhere to its right in the di- 
rection of increasing x (and hence decreasing pressure, 
p, with p — > — oo as x — > +oo). The literature takes p 
periodic in y with period w so the x-values of H at y=0 
and y=w agree. The edges are impenetrable with no sep- 
aration of fluid from them, so that the velocity of fluid 
along them is purely in the x-direction, and so, y=0 and 
y~w are streamlines for the flow. Periodicity of w, in 
synchronizing the flow at both edges is clearly a higher 
symmetry than required. It is readily seen that this im- 
plies center-line symmetry, with v also in the x-direction 
along y=w/2. This, however, is just the symmetry of 
the observed finger, and so the symmetry in which to 
seek a theoretical result. This then is the choice of the 
literature. 

Implicit in this choice of symmetry is, of course, the 
underlying fact that the flow is along the x-direction at 
each wall. This has a striking consequence not noticed in 
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the previous literature. The conformal machinery utilizes 
an invertible analytic mapping £ = h[z) from the fluid's 
spatial plane z = x + iy to a new plane £ = —p + is with p 
pressure, and s the "stream function" , where the curves 
of constant s are always orthogonal to those of constant 
pressure. The thin limit of the viscous fluid equations 
is Darcy's law: The velocity is the negative gradient of 
pressure. Thus, fluid flows with a velocity at a point, at 
each moment of time, in the direction of a curve of con- 
stant s through that point at that time. (These curves 
all evolve in time as then does h.) In particular, with 
velocity in the x-direction along the walls, this is to say 
real in the complex description, since Darcy's law reads 
v x —iv y = h'(z). This implies that at each instant of time 
the part of the strip between y—0 and y=w all to the right 
of H is mapped by h into a strip between s=0 and s=w 
all to the right of p=0 by the obvious choice of the origin 
of pressure (only differences matter). But then, the real 
z-axis is mapped by h to the real £-axis. This, in conse- 
quence of analytic continuation, implies that h is Schwarz 
reflection symmetric, h(z) = h(z). (This, together with 
an identical statement for the upper edge implies period- 
icity over 2w, the true symmetry of the flow. Increasing 
the symmetry to w periodicity is equivalent to an extra 
reflection symmetry through y=w/2.) 

The equation of motion for the interface H is that it 
flow into itself under its velocity at each instant as de- 
termined by Darcy's law. Expressed in the £-plane, this 
becomes a relationship along H (which is just Re£=0) 
of both the C and time derivatives of / = hr x and their 
complex conjugates. With / reflection symmetric, these 
conjugates are then just the functions themselves evalu- 
ated at the conjugate of £. But, along H, the conjugate 
of £ is simply — £. Thus, along H the equation of motion 
is an analytic partial differential equation in £ and time. 
But then it analytically continues as that same equation 
as a field equation throughout fluid and its continuation. 
p3|,p| That is, there is a new symmetry of the equations 
of motion, parity in £. The symmetry of s — ■> — s is the 
periodicity of the literature. 

Unnoticed in the literature, but implicitly present, is 
the symmetry of p — * —p. That is, the equations of mo- 
tion express a relationship between positive and negative 
values of p. This is extraordinary, since it expresses some 
connection between the shape of H, determined by sin- 
gularities in / at (positive) high pressures to behaviors 
of the fluid at (negative) low pressures, i.e. to its efflux. 
This paper explores that connection, and discovers it be 
replete with consequences. 

In particular, for a propagating finger, A=l/2 is 
uniquely determined for the finite geometry with pres- 
sure fixed on each interface just within the viscous fluid, 
and so with zero surface tension on both H and L. This 
immediately implies that the limit I — > oo is singular, and 
S-T pattern selection is unambiguously resolved in finite 
geometry. 

The class of S-T solutions for arbitrary A in the infinite 
problem is particularly simple. When we consider how 



they transport fluid particles arbitrarily far downstream, 
it is easily discovered that only the A=l/2 solution has 
cross-channel curves of particles all at a fixed value of 
pressure all flowing into other curves of fixed pressure at 
later times. With no surface tension, the fluid may be 
opened up to atmosphere along any such curve, which 
then is L. The actual form of the A=l/2 solution came 
from mathematical simplicity in the infinite problem. In 
the finite version it is the unique solution to the ques- 
tion of a translating finger. The other S-T solutions for 
A < 1/2 far downstream (far to the right) correspond 
to solutions for the finite problem with H still without 
surface tension, but with a/v = 1 — 2A on L, the curve 
that can then be opened to atmosphere. This confirms 
the significant connection between the shape of H and 
some behavior (smoothing) of L. Moreover only A < 1/2 
is physically allowed, but requiring enormous values of 
a to substantially decrease A from 1/2. For viscosities 
and high velocities of the experimental literature, this 
amounts to no more than a several percent reduction in 
A. Reference H found just such an order of magnitude 
reduction from 1/2, which there was interpreted as the 
absence of any genuine selection, and speculated that it 
might be in consequence of films of fluid left behind on 
the glass plates behind the advancing finger. Moreover, 
it violates the singular limit theory of a — > + , which pre- 
cisely determines just 1/2. Finiteness can possibly solve 
this puzzle in the literature. 

Next, with a=0 on H for the infinite literature, all the 
S-T solutions are linearly unstable. The a — > + theory 
stabilizes its A=l/2 solution. The finite problem without 
surface tension, significantly, has its A=l/2 solution sta- 
ble in a somewhat anomalous sense. When we consider 
fluctuations about a given solution, we presume that they 
can be imposed with an arbitrarily small expenditure of 
power. In the finite problem the same instabilities, that 
in the <r=0 infinite geometry grew exponentially, require 
power exactly proportional to their size in order to be 
sustained. If there is no source to provide this power, 
then the solution is fully stable in finite geometry. This 
means that if the pump driving the flow is prevented, say 
by some feedback control, from engaging in rapid varia- 
tions in the power requested from it, then the solution is 
stabilized. This is so different from the infinite geometry, 
that we again realize / — * oo is singular. The I — > oo 
limit of the a — > + theory has never been worked out, 
let alone questioned. It clearly should be. 

The idea that the treatment of an arbitrarily distant 
efflux can modify the shape of an object has never been 
experimentally checked. It is worth bearing in mind that 
aggregation arises in diffusive processes for which again 
there is an incompressible fluid with identical dynamics - 
namely the probability flux of random walkers. There is 
enough sufficiently delicate in the considerations to fol- 
low in this paper, that one can only go so far as to spec- 
ulate that quite to the contrary of the accepted wisdom, 
patterns might more generally depend on details remote 
from the selected shape itself. 
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This paper is laid out as follows: 

In Section 2. we erect the conformal field equations 
suitable to two free interfaces, and notice that they are 
displacement invariant under a time (p, the integrated 
flux up to time t. In order to determine exact solu- 
tions to perturb about, we ask the natural questions: 
Are there purely translating solutions, and are there so- 
lutions which automatically satisfy the equations of both 
interfaces, simply in consequence of this translation in- 
variance? 

In Section 3. we establish that the only such solutions 
turn out to be the 1/2 width Saffman- Taylor fingers, here 
not guessed at, but uniquely determined. We determine 
that Saffman- Taylor fingers of a narrower width corre- 
spond to gathering and smoothing the efflux in a precise 
way, and no matter how distant that efflux is from the 
driven interface. 

Section 4. develops the machinery to explicitly solve 
the finite problem in perturbation theory. A felicitous 
symmetry of the equations of both interfaces allows us 
to fully integrate both equations into a purely algebraic 
problem. Although it is known that the infinite case is 
fully integrable, that is not known in the finite geometry 
we discuss. This symmetry allows one to speculate that 
this problem too is completely integrable. A result of this 
effort is that all unstable modes, and only these, deter- 
mine (unstably growing) external impedances, and this 
one function, just of time, uniquely encodes the entire 
unstable spatial flow, thus allowing a pump driving the 
flow to "hear" the precise shape of the interface simply 
by detecting the power requested from it, and so, render- 
ing it capable of completely stabilizing the flow. No such 
thing is true in the infinite case, and so we realize that 
just how the flow is terminated matters, no matter how 
long the body of fluid. 

In Section 5. we realize that while formal pcrturbative 
pole-dynamical solutions |9j-|l^1 that apparently produce 
different width fingers exist, these solutions are purely 
formal, and can't correspond to any actual approximate 
solutions in the finite problem nearby to the 1/2 finger. 
This fully establishes isolated pattern selection. 

The Discussion Section 6. more iconically presents 
what this paper has accomplished after the reader has 
been armed by the machinery presented. 

The paper ends in an Appendix that establishes that 
the pure S-T 1/2 finger is unique, and so totally isolated 
within pole-dynamics for the finitely terminated prob- 
lem. It has been relegated to an appendix in order not to 
interrupt the main line of exposition. This Section con- 
tains some modest new results for pole-dynamics, and 
is a self-contained extension to paper §. That paper 
can also assist the reader in gaining familiarity with the 
workings of our reflection symmetric field equations. In- 
deed, some of the material of this Section was reflected 
into that paper. 

As a final point of information, in each numbered Sec- 
tion, e.g. 2., equation numbers start again from 0, and 
are referred to within just that same Section without a 



Section number prefix. References to equations in Sec- 
tions external to a given Section are then prefixed by the 
number of the Section they appear in. 



2. THE GENERAL SETUP AND EQUATIONS OF 
MOTION 

Taking z = x+iy, with x increasing down-stream along 
a channel of cross- width y G [0, it] , the fluid obeys Darcy's 
law for the velocity v and pressure p 



v = —dp 



(1) 



The flow is incompressible so that d ■ v = = Ap, and so 
p is harmonic, and the real part of an analytic function 
h, naturally constructed through its derivative: 



ti 



iv y -d x p + id y 



V(t) V(t) 

-p + is 



V(t) 



h — £ + irj = 



v(t)) 



(2) 



which satisfies Cauchy-Riemann for Ap = 0, and s is the 
harmonic conjugate of —p. V(t) here is a function of time 
alone. Consider the change in h along an arbitrary curve 
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8h 



1 i 
h'dz = ——5p+— I v n d£ 



(3) 



Since by Eq.(|) 6( = -Sp/V(t), 

Sri = W)^l Vndi 



and so 



Srj = 7T for V(t) = - [ 

71" Jy 



At 



(4) 



Thus with V(t) the mean velocity of all fluid (the con- 
served flux divided by cross- width 7r), the original spa- 
tial channel is analytically mapped to an identical one 
in £. (Although flux is conserved throughout the incom- 
pressible fluid, for a finite system it will generally vary in 
time.) 

Assuming the fluid not to stagnate other than in front 
of added boundaries, d / 0, A' / 0, and so h is invertible, 
conformally mapping one [0,7r] strip onto another. We 
denote its inverse as /: 



f = h-\f 



V 



at z = /(C) 



(5) 



The virtue of / is that all boundaries are known in 
Q: The problem we consider has impenetrable walls at 
rj = 0, 7r (along which v is real) and two free interfaces, 
one at p = 0, i.e., £ = and the other at p — p a tm = ~Pg, 
or at 
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(6) 



p g > is the gauge pressure across the body of fluid, 
and £ 9 (the total impedance) positive and generally time- 
dependent, because either or both of p g and V are. (p g is 
time dependent if we experimentally modulate the driv- 
ing pressure.) Thus, in C, the fluid is just a rectangle 
extending from £- = up to a variable right hand edge. 

To track a fluid particle through use of the one- 
parameter (time) family of maps, /, tag a particle at 
t = by its ^-coordinate Co : 



Then 



*(Co,*) = /(C(G>, *),*); C(Co,o) = Co 



zt = f'Ct + ft = V/f\ or 

v = \f\ 2 Ct + hf 



(7) 



(8) 



Once / and V are known, Eq. (|8|) is an ordinary differ- 
ential equation to be solved for C(Co , given the initial 
data of Eq.(0). 

To determine / requires the physical input that an 
interface simply transports under its velocity. On an in- 
terface ReC = £(t), with £ = on the left interface and 
i = €g(t) on the right. That is, if ReCo = £(Q) then 
ReC(Co,£) = £,(t) for a fluid particle on the interface. 
That is, 

ReC* = i(t) on ReC = , (9) 
so that (||) is 

V(t) = Re(/'(/'i + f t )) on ReC = | , (10) 

on each interface. The form of (10) recommends the fol- 
lowing definition: 

= + (11) 

immediately yielding 

V(t) = Re(/|/| t ) on ReC = (10') 

on each interface ReC = 

The generally unknown-in-advance V(t) complicates 
(10'). Delaying its determination to a final step, it is 
overwhelmingly expedient to now eliminate it through a 
redefinition of time: 



<P = V, 



(12) 



so that tp-time is integrated flux to time t, that is, the 
volume of fluid moved, or the energy expended by exter- 
nal sources in the case p g = const. But then dt — Vd v 
and (10') is simply 



Finally since / takes ImC 
reflection symmetric, 



But ReC = 



/(£*) = /(£*) . 
C = -£ and (10") reads 



to Imz = 0, / is Schwarz- 
(13) 



2 = f'A-QhJO + f'AOfi.J-0 • 



(14) 



Although (jTJ) is exactly (10") on ReC = 0, ([jj) is in 
fact correct over all of /'s analytic continuation, since 
(14) is an equation just of the analytic variable C and 
(14) and (10") agree all along the curve of an interface. 

For our problem of two free interfaces, let us call /^_ 
simply /, and fa = g: 

S(C,¥>) = /(C + ■ (15) 

Then the equations of motion ( |l4"| ) are simply 

2 = f\-0U{0 + f'(QU(-Q , (14') 



ff , (-o^(o+s / (o^(-o 



(14") 



1 = Re(/|4 v ) on ReC = 



(10") 



That is, / and its £g translate g are both solutions to 
(14'). Although (14') is strange in its C-dependence, it is 
simply autonomous (translation-invariant) in ip. That is, 

/(C, <f) a solution to (14') ^> f((, f + ifo) also a solution. 

(16) 

This is significant in that purely translating solutions 
are allowed and a path is also opened for both f and g to 
obey (14') in a way independent of the matching-up of 
inner details, merely reliant on translation invariance. 

Pause to notice that reflection symmetry for / has 
determined a new symmetry, parity, C ~ * f° r t ne 
equations of motion (14'). This entails not only a cross- 
channel symmetry, but also a relation of p to —p. But 
then, far upstream singularities which determine the 
shape of the driven interface become related to far down- 
stream properties, which is to say how the flow is termi- 
nated. This observation is the foundation for the work 
that follows. 

Not only is / reflection symmetric about the £-axis, 
but also about ImC = T since the upper horizontal wall 
there is also impenetrable. Coupled with f's reflection 
symmetry, this is simply the periodicity of /: 

/(C + 2t™) = /(C) + 2th. 

Reflection symmetry itself means that both the chan- 
nel and its C-image may be mentally reflected through the 
x and £-axes respectively, so that, in light of the above, 
both views are periodic in y G [— tt, 7r], although only the 
upper half is physical. However, the experimental flow 
actually is very close to symmetric about its center line. 
We now capitalize on this by taking the physical channel 
to be [— 7r,7r], so that reflection symmetry now enforces 
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actual symmetry, and later, if we so choose, study the 
stability of this symmetry. 

That is, our channel is now of width 2-k and Eq.(||) is 
for the sequel modified to 

V(t) = J v n dl . 



3. PATTERN SELECTION AND PHYSICAL 
FINITENESS 

Since the equations for the interfaces are ^-translation 
invariant, within an available infinite channel we can seek 
solutions for which the driven (j> = 0) interface retains 
its shape, but merely translates in ip-time. By analytic 
continuation, such an / must be of form 



/ = /%)+u(C). 

Then, by (2.14') 

1 _ u'(Q+u'(-Q 

n<p) 2 



(i) 



(2) 



for A some real constant. Neglecting any branching con- 
siderations and taking u analytic at Q = 0, (||) integrates 
to 



u(C)-u(-C) 



AC 



(3) 



for the antisymmetric part of u. On £ = Q , with u 
reflection symmetric 

Im u{irf) = 2/(77) = Ar/ 

so that the driven interface occupies a fraction A of the 
channel width. Thus A < 1. With A < 1 the interface 
must stretch off to minus infinity. By "finite" we shall 
then request that at least the right interface (the "efflux" ) 
be of bounded re-extent. These solutions with A < 1 are 
called fingers in consequence of their long asymptotically 
horizontal sides at ±i7rA. (We will determine truly finite 
configurations at the end of this section.) 

By §), 1) 



/ 



"(0 



with 



U (C) = AC + £(C); E(-0 = E(C). 



(4) 



(5) 



Clearly E must have a discontinuity of 2m(l — A) 
throughout fluid to meet channel wall data, so that E 
must contain logarithms. Otherwise neither A nor E is 
determined. 
Next 



= ! + «(£ + &>)) 



must again satisfy (2.14'): 

2 = (u' + + u'_)/\ + 2£'*4«'_ with u'± = u'(±( + £). 
Multiplying by 2X 2 /u' + u'_ and rearranging, 

F(t+0F(£-Q = l + 4\ 2 i'(tp) (7) 

with 

F(C) = 2X/u'(0 - 1. (8) 
Substituting (^) into the second of (||), 

F(C)F(-C) = 1 (9) 

so that 

F(0)=±1. (10) 
Setting C = i in (|^), at which it must be defined, 

F(i + QF(i-Q=±F(2g) (11) 

and 

1 + 4A 2 C' = ±F(2C). (12) 

There are two possibilities: C varies with (p or is just a 
constant. If £' = then (Q) and (||) become 

F(£ + 2g) = F(Q 

and hence u' is doubly periodic with real period 2£ and 
imaginary period 2ni. But then, by (||) E'(() is an odd 
such doubly periodic function. Not to be just constant 
(and hence zero), E' must have singularities in each cell. 
They cannot lie in £ 6 (0,C) where / must be analytic. 
But then, they can't lie in (£, 2£) either, since then, by 
periodicity they lie in (— £, 0) and then by oddness, in 
(0, C) again. Thus singularities must appear on cither 
C = or C = C- With A < 1, there must be simple poles 
to pick up the requisite 2ni(l — A) discontinuity. With no 
net integral around a cell, the sums of the residues must 
vanish, and so there must be at least two poles per cell. 
Finally, the cut in E must run the length of fluid, and 
so there must be poles on both £ — and C = C- But 
then, the right interface must also be unbounded, which 
we have rejected by finitencss. That is, we have 



A = 1, f = tp + ( 



(13) 



or £' f 0. 

If £' ^ 0, the only continuous solutions to ( fj"l| ) with £ 
and C independently varying is exponential: 



(6) 



F(() = ±e-" c 



(14) 



G 



with n an integer to meet 2iri periodicity. But then, 



and so (2.14') for just the driven interface is satisfied by 



2A 



1 ± e""C ' 



O \ 1 

u = — ln(e" c ± 1) = 2A(C + - ln(l ± e~ nC )). (15) 



Thus, n > 1 (fluid to the right of £) and X — 1/2, so 
that 



/ = 2^ + C+-ln(l±e-" c ) 
n 



(16) 



and wall boundary geometry is satisfied, since the coef- 
ficient of £ is 1. n > 1 is simply n parallel fingers each 
of width l/2n, and for simplicity, and the one finger case 
we care about from experiment, 



/ = 2(^ + C + ln(l + e- c ) 



(17) 



where we have chosen the + to have the nose of the fin- 
ger at channel center. This is precisely the 1/2 width 
Saffman- Taylor finger. However, both A = 1/2 and the 
precise form of E are here totally and uniquely deter- 
mined. That is, with the efflux physically treated (the 
right interface at atmospheric pressure) pattern selection 
is complete. (The great bulk of this paper consists in 
comprehending its stability.) 

Returning to (|l^) to determine £, 



I' = -l + e 



or. 



or, with n = 1, 



2n£ _ 



= l + fce 



-2ntp 



=2? _ 



1 



ke~ 2v 



(18) 



where k determines, at reference time ip = 0, the length 
io of fluid from the center of the driven interface at In 2 
to the center of the right interface at ln(e^ < -°- > + 1), or 



k = 4e 2L °(l 



-L 



(19) 



Before further considering same details of this solution, 
and rendering the solution truly finite, let us consider 
what we have accomplished, and why efflux treatment 
determines selection. We proceed in two steps. First, we 
consider a class of "solutions" containing ( p7| ) in order to 
see just how sharply A = 1/2 is selected, and how it is 
contingent upon the treatment of efflux. In the bulk of 
this paper that follows, we consider if ([!]), exact trans- 
lation, isn't simply too precise and perhaps not quite 
physical a request. (We shall vindicate it.) 

By @ and (|), 

E = ln2coshC/2, 



/ = p/X + A( + 2(1 - X)E 



or. 



/ = ¥>/A + C + 2(l-A)ln(l + e-t), 



(20) 



for any A G (0, 1). These are the Saffman- Taylor fingers 
for width A. As we already know, there is no £ g for any 
A / 1/2 for which g will also satisfy (2.14'). Let us see 
why. 

To track fluid particles under (|2^), we determine 
C(Co,¥>) by (2.8): 



(e' + |l + (2A - l)e^ c | 

= (- - 1) (e" 2€ - 1 + 2ie-«sin?7) 
A 



(21) 



It is immediately clear that the equation decouples in £ 
if and only if A = 1/2, in which case 



-2? 



1 



(just (jig), of course) and fluid particles all with the same 
value of £ at time ip, but with all different rj values, flow 
again at later times into all with the same value of £. This 
is so no matter how far downstream we look. And just 
so, for A 1/2, no matter how far downstream, a curve 
of fluid particles all at the same £ at time ip will fail to be 
so at any later time. That is, by opening the channel to 
the atmosphere all along the curve £ = £((p) for A = 1/2, 
the flow will continue to remain at atmospheric pressure. 
However for A =/= 1/2, this flow may never, no matter how 
far downstream, be able to be opened to atmosphere. 

But with / differing from free flow by 0(e~^), which 
for large £ is astronomically smaller than any physically 
measurable entity, perhaps our question of A = 1/2 is 
purely mathematical, and not physical. Let us now see 
that it is highly physical, with measurable consequences, 
indeed consequences that conceivably have already been 
experimentally measured, and in the literature. 

By© 



(If) 



smb. £ 
sinry 



and so, 



tanh — tan — = tan — 
2 2 2 

for that curve £(</?) which at £ 
(E2fi has the series development 



(22) 



-oo ends at n = 



77 = + -e _< sin n0. 
^ n 



(23) 



We proceed now to solve 
taining terms just of 0{e~ 



U) just far downstream, re- 
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+ 2(2A - l)e"« cos??) - -(i - 1). (24) 

A 



By (23) to 0(1), 



COS T) ~ COS ( 



and integrating, 



e-2(2A-l)e-*oos0~£ o (0)-(^-%MMKo » 1 
and so, 

£~£( v , ) 0)+2(2A-l)e^cos0, (25) 
and by @, 

77 - 6» + 2e"£ sin (9. (26) 

Next, to 0(e"«), 

z = / ~ £ + 2(2A - l)e- c + ip/X 
and so, substituting (p5|), (|26|). 

2-eo + ^ + «^ + 2Ae-« +ie , (27) 

and at each time ip, 9 parameterizes a curve of fluid parti- 
cles spanning the channel. For £o large enough, at tp = 
these are just a flat interface by selecting those fluid par- 
ticles with £o(0) = £o independent of 9. 

Let us now compute the curvature of z(i9). This is 
just, with the center of the curvature to the left, 



1 z" 

K = W\ \~' 



With z' = 1 + 2Ae- ?+ie , to 0(e"«), 
k = 2Ae~^ cos 9 

and so, (f25f) reads 



£0) ~ € + (\ ~ 2)«. 



(28) 



(29) 



Consider now treating the efflux, namely endowing the 
right efflux interface with a surface tension a, and main- 
taining the pressure against this interface constant along 
it. Then, 

p((p) = p — ok 
and dividing by —V(<p), 



(30) 



But, this is precisely ( (29| ) with 
A ^ V 



and this result no matter how asymptotically far down- 
stream we open the flow to atmosphere. Now, /' of j2fj| ) 
at ( = (the nose of the finger) is A, and so, the mean V 
is related to the finger's speed v by V — Xv. Then (|30| ) 
is 



1 - 2A 



v 



{2n) 2 B ~ AOB 



(31) 



with the dimensionless B that of the literature 

, 2 



i,-i = i2^ ry 

a V o / 



with /j viscosity, w/b the ratio of channel width to the 
gap between the upper and lower plates, and 2tt comes 
from our natural choice of channel width, rather than 1. 

We now see that with a > 0, we can only produce fin- 
gers of width smaller than 1/2 by physical smoothing and 
"gathering" means, and that the surface tension must be 
physically enormous to significantly decrease A from 1/2. 
For the physical choices of a at experimental conditions 
when one observes a/v so small as to be close to zero 
surface tension, one has B ~ 10~ 3 . Had this been our 
configuration, with no surface tension on the driven in- 
terface, but this value for the efflux interface, we would 
produce a A ~ .48. Because we now realize that how 
the efflux is treated does matter, and since we have no 
idea just how it was accomplished in the published lit- 
erature (because the experimentalists presumed it didn't 
matter), all we can say is that we have explained both 
qualitatively and to within order of magnitude a puzzle 
in the literature || . 

Just as importantly, we have determined that the infi- 
nite channel-full of fluid is physically an incomplete spec- 
ification of an experiment. In never physically treating 
efflux, the pole at Re £ — > +oo of ( |20| ) promiscuously 
allows the physical completion of an experiment in any 
way whatsoever. For a given physical method of collect- 
ing efflux, the fluid can be opened to air at any large, 
but nevertheless, finite length down the channel. This 
limit can then be taken to infinity, but for example, in 
our configuration with a second interface, leaving intact 
the selection of (|3l|). That is, the pole at infinity is ac- 
tually, and physically illegally, within fluid, and so the 
usual theory of the infinite channel is too compliant to 
diverse stresses. 

Since finiteness matters, it is useful to determine a fi- 
nite equivalent to ([l]). We can utilize tp-invariance more 
deeply by requesting that the second interface is none 
other than the first, but delayed in time by ipo: 

/(C + i(<p, <po), V) = f(C, V-Vo)+ /%o). (32) 

By this we mean first that a constant, 0(tpo) is allowed 
since (2.14') entails just derivatives. Clearly, the amount 
of delay is related to the length of the body of fluid, so 
that the impedance £ at the right interface depends upon 
it as well. We seek solutions for any chosen length of fluid, 
and so for arbitrary ipo . Finally, we seek a shape of the 
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driven interface that, in the spirit of (P is independent 
of just what length we contemplate so that / itself is 
independent of ipo . As the length of the body of fluid 
goes to zero, as with | of (|l^) with k — > 0, we request 



£(p,0) = =► 0(0) =0. 

Let us now determine the solutions to (£52| 
Set ipo = <p in (p2|) and define 



then. 



(33) 



S(<p,<p) = M<p) ; /(C,o) = M (c). 
/(C + C.(p),v) = /9(v) + «(C) 
/(C,^)=/3(v)+«(C-C.(v))- 



(34) 



Substituting in (J32|) , with 

j(<p,<Po) = €(v,<po) - C*0) + C*(<£ - Vo), 

obtain 

u(C + 7(V. ¥*>)) - "(C) = /% - ¥>o) - /%) + /%o)- 

Differentiating on £, should 7 vary with y>o then it' = 
const, and u = £. Otherwise, with 7 constant, setting 
ipo = 0, by ( |33| ) 7 = 0, in which case /?(y> — <^ ) = (3(<p) — 
fl(ipa) and so /3((p) — y/A for some constant A. Thus, 



and 



/(£?>) = f + u(C-C.foO) 



|((/J, V? ) = C*0) -C*(V-^o)- 



(35) 



(36) 



We now ask that ( |35| ) satisfies (2.14'). But this is pre- 
cisely the discussion for g of (@) upon replacing £(92) by 
— C*(^)) but without (H) so that F(0) = a is unknown. 
Certainly ^ 0, since otherwise f = by (|3^). So, by 

(ID, 

«(£ - £„) = 2A(C - C* + - ln(l + ae-"K-W)) 
n 

with ae _Tl ^ _l >*^ = ie^'^-C) and C < 0, we must again 
have A = 1/2 for channel geometry. ((IJ) is now 



C = 1 - a 2 e 



2„2<, 



or, 



C' = l 



together with 



/ = 2^+-ln(e"( c -^±l). 



(37) 



But this is just the simplest pole-dynamical 1/2 finger, 
converging to (16) as <p — > +00. Indeed, ( fl6| ) is the limit 
of these finite solutions as the length ipo — > +00 and 
we take <p — > +00, but leaving a finite length of fluid 
downstream. 

Thus, with n = 1 and the + sign, and in the above limit 
of a truly finite problem, we have precisely two solutions 
related to translation: 



A = 1, / = ip + £ i = La = const, e 2 ^ = k 
and 

A = 1/2, / = 2tp + ln(e c + 1), e 2 « = 1 + ke 
For a fully finite version, we have 
A = 1/2, / = 2tp + ln(e c ^ + 1), e" 2 ^ = 1 + e 



-2tp 



lip 



i = C(V>) - C(f - fo), e 



2€ _ 



1 + ke~ 2ip 
1 + e~ 2 v ' 



(38) 



with 



(po = L = A/2ir, A = area of fluid, k 



o 2 VO 



4. REFLECTION SYMMETRIC 
PERTURBATION THEORY FOR TWO 
INTERFACES 

Having asked, in the light of experiment, for purely 
translating solutions, we must now answer the question 
if the fluid in this theory of no surface tension actually 
assumes the selected solution, or instead, even when ini- 
tially set, deeply unstably runs away from them. There 
is a related question to this. For Lq 3> 1 in (3jl9|) and 
p « Lq, by (3|D 

f ~ L - ip + In 2 

and £ together with the length of fluid downstream from 
the driven interface are very large. By the time there 
has been a net flow of fluid ~ Lq, £ ~ 1 and the right 
interface has begun to bend into a 1/2 finger parallel 
to the left interface, with the distance between the two 
interfaces rapidly turning into a thin film for larger values 
of ip. It is intuitive that a long body of fluid downstream 
can help stabilize the driven interface, but certainly not 
just a film of fluid. That is, one can expect qualitatively 
different behaviors with « Lo and ip ~ Lq . Let us 
call this former regime the "finger regime" . This raises a 
new question. Is it possible that during the finger regime 



we have solutions exponentially close (i.e. 



b «) to 



purely translating, but then rapidly, for tp ~ Lq becoming 
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highly dynamical. Then, from the physical viewpoint 
of an experiment, we observe a "perfectly" translating 
finger for the duration of observation. This makes (3{l]) 
again too sharp a question, and the theory might well 
allow for /'s similar to those of (3|2(]) for arbitrary A, 
but decorated by dynamical corrections of 0(e~ L ") over 
the finger regime. 

In a way, such a possibility is curious. It implies that 
for (3.pp|) a powerful surface tension is requested on the 
efflux, while (3.|o|) decorated by 0(e~ L °) corrections pro- 
duces exponentially close motion throughout the great 
bulk of the fluid, while requiring no surface tension at 
all. This is surprising since by (3.^5|), for any A, 

£ ~ L Q - (I - IV, 

so that macroscopically different fluxes are elicited with- 
out physical agency. 

Regrettably, beyond our A = 1/2 exact solutions, we 
know of no way to construct any other exact solutions 
with two interfaces. For example, it is not difficult to 
show that the pole dynamics of the literature is exhausted 
by just (3.|38|). (Since the details of this excersize are tan- 
gential to our main thrust here, it has been relegated to 
an appendix, Section 7.) This leaves us in the unfortu- 
nate position of a merely perturbative treatment, which 
we now erect. 

Perturbing a solution to (2.14') determines a very spe- 
cial change of variables which renders the calculations 
feasible. Donate by / a known solution to (2.14'), call its 
perturbed value / p , the strength of perturbation e, and 
/ the perturbation: 



dip 



f P = f + ef + 



(1) 



We shall denote by "transpose" the symmetry operation 
of (2.14'): 



so that (2.14') for / is 



/'% + (/'%)' - 2. 



(2) 



(3) 



Entering (g) into (2.14') leads to 0(e): 



or, 



/'% + f v f + transpose = 0, 



/'% + /£/' = a(C,^); a* = -a. (4) 



Integrating (Q) on curves in a solution surface, with s a 
parameter, 



g = /'(-c- dC 

By the first two, 



5 = /'(-C^),g = /,(-C^),|=a(C^). (5) 



o = ~ f'(-C<p)-n = -rfi-CM- (6) 



ds 



ds ds" 



That is, we should take a function purely of /* as one new 
variable, and take a function of (£, tp) independent of /* 
as ,s. Then, the partial derivative of / with respect to this 
second new variable is, by the third o.d.e., antisymmetric 
in £. This optimally selects the transpose of the first 
variable as the second variable. To ward off exponentials 
when we consider an / such as (3.0), we are led to define 
the following two new variables in place of (£, ip): 



v 



e f(-C,<p) ) £ = e /(C.¥0. 



We then have by (|2|) for any function u, 



(7) 



(8) 



the reason we have denoted the symmetry of (2.14') as 
transpose. The last of @ with s = In £ is 

since the mixed derivative of an antisymmetric function 
is also antisymmetric, and so 



/ = vd v a{v,£,). 



(9) 



(The integration constant, a function of v is just an 
a(y, £) = u(v) — u(£).) That is, the solution to first order 
perturbation, (^) are In v derivatives of an arbitrary an- 
tisymmetric "potential" a, and first order theory for one 
interface has been fully integrated, despite the non-local 
(in C) form of the p.d.e. (2.14'). 

Let us check that our new variables v, £ are sensibly 
defined. Now, £ is complex while ip is real. Both v and 
£ are complex. If they are formally independent, then 
upon inverting they determine a complex extension of (p. 
But they are independent, precisely in virtue of (2.14'): 



<9(ln£, In v) 

d(C,p) 

and so, by (2.14') 

<9(ln£, \nis 



/ ftp 



= 2. 



Inverting fllO|), 



5(ln£,hH 2 \, /" /' 



1 ( fin fif 



(10) 



(11) 



(12) 



Paying attention to (3.|T^), rather than ip, e~ 2ip will 
constantly appear, and so we choose for later a new time 
variable 



A 



-lip 



(13) 



for which we have by the second row of (|l2|), 
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vd v \ = —A/' and £<%A = —A/" 



(14) 



which will shortly become important. 

Before proceeding it is useful to rewrite (2.14') in (v, £) 
variables. By @, 

3< = /'^ - /'V^ , 9 V = /^c\ + . (15) 

Substituting in (2.14') for /, then yields 

^ = hf\ - fuft- (is) 

Clearly / = In £ obeys (|l^) , while i In v is not reflection 
symmetric. It is now trivial to obtain the perturbation 
results, which to second order are 



and so, 



and 



/ p = ln£ + e/ + e 2 F- 



^ + 7/1 = 



(17) 



(18) 



(19) 



( |l8| ) is, of course, the observations below (||), and its 
general solution just ([)]). It is also easy, by successive 
integrations by parts and the use of ( jig ) to integrate 
@ fully into 



all with p > unstable. Here we can also take for the v 
term of a J ^ In v producing 

/ = lnz/ = 2<^ + ln(l + e- f ) 

and so perturb to the nearby Saffman- Taylor solutions 
(3.po|) with A = (1 — e)/2. Indeed, any function w(^) 
analytic at Re£ — » +oo provides a legal perturbation. 

Let us now go on to the second interface. Not only 
does / satisfy (2.14'), but so too does 



g(C,<p) = f(( + 

Accordingly, we define 



(21) 
(22) 

so that, by @, 

HC, <P) = £(C + I, ¥>) and P(C, ¥») = ?(-C + I, ¥»). (23) 
It now follows that 



g p = \n( + eg + e 2 G 



and so by (||) 



g = vdua g {v,l), 



and by (EC 



G = . 9 |%g + i>a iJ (-«?« 7 t + A 9 ) 



(24) 



(25) 



(26) 



(20) 



with another antisymmetric potential. 

So far we have just discussed one interface. This suf- 
fices for the infinite channel. With the flat solution, 
(3.0), 



and so 



v v CP 

— => f = v p 

P P 



-pC+pv 



with p the positive integers is a full basis of perturbations 
for the one interface at Re£ = 0, and all are exponentially 
unstable. 

For the 1/2 finger solution of (3.|l7j) (and almost iden- 
tically for (3.§)), with A of ph, 



Z = 



l 



i 



and again 



v p CP 

a= 5- f = v p = (l 

P P 



Cy e 2pv 



with a g and A g antisymmetric in (£,£). 

But f p and g p are related by a C translation, by (2jl5|) 



(2© 



Here we must pause. Should Zg be just £ that relates / 
to g, or should it be modified, basically in proportion to 
the strength of the perturbation /? Should Z be modified 
to 



Z g = £ + + 



(27) 



(each term a function of ip), then certainly these fluc- 
tuations can not be freely imposed, since under a fixed 
pressure difference, p g , changing £ changes the flux of the 
flow, and so costs power of an amount proportional to the 
size of the fluctuation, and makes it possible for the exter- 
nal agencies driving the flow to control the fluctuations. 
Should ^, etc., however, vanish then the corresponding 
fluctuation is free and uncontrollable. 

To be general, substitute ( |27| ) into (2.|l5|) and expand 
in e to relate g and G to / and F: 



(28) 



1 



G(C, (p) = (F + /'$ + *3 C (/ + -m)(C + Z, <p) (29) 
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where = since $ is a function just of ip. 
Start with @: 

<?(C - I, *0 = / + /'*■ 

Let us define for convenience 

* = A^'(A). 



Then, by (|14|) 

/'* - -vd v \i>'{\) = -vd v ip(\). 
But then, by (g), 

^(a - ^(A)) = ,g(C - |, </?). 



(30) 
(31) 
(32) 



This scheme succeeds because of an extraordinary prop- 
erty of the (y, £) variables. 
By m 



6* =£(C-?»=£(C, and 



(33) 



So, by (|llj), valid just as well for (£, £), 
<9(ln £*,ln^») 



0(C,v) 

But then, dividing by (|1 
d (In £*,lnz/* 





9(ln£, ln£) 




l-l' 








1 



= 2. 



d(ln£,lni/) 



= 1 = 



<9(ln£, ln^* 



<9(ln£, In v) 



dim 



d In v 



and so, 



^=£(2£-C,¥>) = f/7(0- 



(34) 



The remarkable result is that 7 is a function purely of £. 

Let us elaborate on The simplicity of its deduc- 

tion reflects the curious fact that the equations of motion 
are just the constant Jacobian of (/(£, <p), /(— C, y)) with 
respect to (Cj¥>)- ^ i s straightforward to directly deduce 
it by the joint integration of (2.14') for both f((,tp) and 

ff(C,¥>) = /(C + ■ 

Writing 



7(e) = e r(inc) =e r(/( C , v )) ; 
assuming (|34|) we have 

/(-C,v)-/(^-c.v) = r(/(c, ¥ >)) 



(35) 
(36) 



where again, it is remarkable that T is a function purely 
of/(£v). 

To directly deduce this result, write down (2.14') for 
g = /(C + £s) an d translate the argument £ — > £ — £ 9 . 
Subtract from this (2.14') for / and so obtain 



-/„(0W(-C)-/(2£ fl -0). 

Dividing by J'/^, the equality of ratios is equivalent to 
@. That is, if just / obeys (2.14') and satisfies (§|), 
then / is a solution for the finite problem, and conversely. 
This will enable us to put severe restrictions on the na- 
ture of solutions near to /, a pure 1/2 finger solution; 
sufficiently severe to prove pattern selection. 
Now, by @, 

and so, 

9((-£, <P) = ^*a 3 ,i(^*,C*) 



7 (0 a " 1 S(0 ,fJ 



^ a s(-77T»0- 
71?) 



But now, by ©, 



o(i/,0-^(A)-K0 = a fl (rgy,0 



(37) 



(38) 



with the integration constant p purely a function of £, 
and we have succeeded in simultaneous integrating both 
equations of motion. 

At this point a few observations need be made. First 
A is independent of £• Writing \{y, £), the value is un- 
changed by sending £ — > — £, and so A is symmetric in ^ 
and £: 

\{v,t) = \{£,v). (39) 
(This is implicit, of course, in (0).) Next, write Q as 

7tf(C,y))= i ( r C ' y) ■ (40) 

e(c-2o i/(2| - c) 



But then 



7(".)=7(-77x)=7(£(2£-0) 
71?) 

and so, together with (|3^) , we have 



7(C) 

But then with A independent of £, 



7(0 



A(^)=A(£ 7 ( — ),— ). 

7(6 7(0 

Substituting v — » ^7(£), 

AMO. = ^(frM f) = A(i/, : A* = A 



(41) 



(42) 



(43) 



and A is another symmetric function of v and £. Substi- 
tuting v — ► 1/^(0 m (la) now yi^ds 



= /'(c)9 ? (/(-o - /(2e ff - o) 



(44) 
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Transposing ([ll]) and adding, by the antisymmetry of a g 
we have 

(45) 

and we have determined what a's in (^|) are allowed for 
the problem with two interfaces. 

Let us see that the decomposition of (Q) is unique, 
so that for a given perturbation /, and hence a given 
function a, both ip(\) and p(£) are uniquely determined. 
(Were this not the case then ^ would lose meaning.) 
With A and p denoting the differences of two different 
decompositions, we have 



-2^(A)=p(e))+p(i/) (ie.a = 0). 
But, with a = 0, 



(46) 



is 



and by 



But then 



or 



a»("/7,0 = -^(A)-p(O. 

3).© 

ff(C - D = - /'*• 

ffp(C-D = /+e/'* = /(C + e*), 



9p = f(C + (£ + e*), ¥>) and / p = f (a = f = 0). (47) 

That is, the ^ of ( p6| ) corresponds to the identical unper- 
turbed motion / throughout the fluid, but we have taken 
a new interface at £ p = £+ e^f. This means that we have 
simply considered a larger finite body of fluid, and so ip of 
( [flf ) is the perturbation of changing the volume of fluid. 
This is not what we consider, since the physical experi- 
ment is predicated on volume preservation. That is, for 
the physical perturbations we care about, which contain 
no p art of p of the form of that which is determined by 
([46|), our decomposition is unique. 

At the end of Section 3, in (3.|l|), (3.0) and (3.§|) 
we recorded the fluid volume through the parameter k, 
which is e 2L ° with 2ttLq the finite area of fluid. In chang- 
ing Lq by e, we change k by 2ek and so the \& in (^) is 



dk ' 



(48) 



Thus we have for the flat solution by (3.|l3|) 

*l = 1, ipL = InA (49) 
and for both the A = 1/2 fingers of (3.19) and (3||) 
k\ 



*L = 



1 + kX 



ij}L = ln(l + k\). 



(50) 



Let us now put in evidence the functions A, 7 and A 
for these same solutions. 
For (3.0), 



A = 

7 = e~ 2 $ = l/k 



(51) 



and A 

Similarly, for (3.0), 



A = 



7 : 

and A 
Finally, for dynamical (3 
A 



1 1 

£ v 

i 

1 1 

? + ; 



(52) 



7 

and A 



k 
1 

v 



k 



1 



k 



(53) 



Now, by ( |46| ) and (|49|), ( |50| ) we can determine the for- 
bidden forms of p that generate pure volume changes. 
For (3. pi) 



p L (0 = 21n£-lnA, 
and for both (3.0) and (3.||), 

p i (0 = -21n(l + fe/0- 



(54) 



(55) 



Although we have put these details in explicit evidence 
for the exact solutions we know of, it is to be stressed that 
the entire machinery of this section is directly applicable 
to the perturbation on any f,g which solve (2.14'). In 
particular, it will apply to the perturbations of our ex- 
plicit solutions, now regarded as new fundamental solu- 
tions, up to the order we consider. 



We now need to determine a(i>,£) from (38), or elim- 
inating a g , from ([45|). With a antisymmetric, we can 
always write 



a(v,Q = u(y,X) - A). 



(56) 

Entering ( |56| ) into ( ^5| ) and then substituting v — > 1^/7 (^) 
and rearranging, we have 

(u(v, A) — u(-^-r,X) - p(£) - ip(X)) 

7(6 

+ K£7(-77v), A) - A) - _ ip(X)) = 0. (57) 



W(0 



7(0' 



This equation is much more transparent than it appears. 
Reverting to invertible (£, A) variables, with 
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u(^A)=u(£(-C,A),A) = /i(C,A), (58) 

a(C,A) = /i(C,A)-/ l (-C,A) (59) 
and (H), by © is 

(/i(C, A) - h(( - 2& A) - P (aC, A)) - ifr(A))+ 

(M2| - C) A) - ft(-£, A) - P (C(2| - C, A)) - tA(A)) = 
or, with the definition of r(£, A): 

fc(C, A) - ft(C - 2l A) = p(e(C, A)) + V(A) + r(C, A), 

(60) 



r(C,A)+r(2£-C,A) = 0. 



(61) 



It is now easy to see that r can be taken to vanish with 
impunity. This follows because we are uninterested in h 
itself, but rather by (|59]), its antisymmetric part. But, 
( |60| ) is linear in h, and so can be decomposed into three 
pieces of h: One piece satisfying (6C ) with r = 0, a second 
with p + ip = with r satisfying (61), and finally a purely 
homogeneous part 



and so with 



VC,A) = /i (C-2£,A) 



a (C,A)= ao (C-2£,A). 



(62) 



As we have already discussed, an antisymmetric doubly 
periodic function must have singularities within physical 
fluid, so that ao can be dropped by simply determining 
the rest of h to have correct analyticity. Returning to the 
second r contributed piece of h, h r , and suppressing the 
A-dependence, 

MO - MC - 2|) = r(C). 



Setting ( 2£ — and using (|6_1 

M-0 - M2£ - = r(0 = MO ~ h r (C - 20, 

or 

a r (C,A) = a r (C-2|,A) 

which then can be totally absorbed into ao, and then 
ao dropped completely for an h with correct analyticity, 
which then simply obeys 

h(C, A) - h(( - 2l A) = p(aC, A)) + iKA). 

For a convenience of geometry (what we refer to men- 
tally as Re£ — > +oo) , it is convenient to redefine 



/i(C,A) -» -ft(-C,A), which leaves (|56j) and (|59|) un- 
changed, and produces, after sending £ — > — £ 



h(C + 2|, A) - MC, A) = p(i/(C, A)) + V(A). 



(63) 



- the one which possesses 
is then the unique solution 



The unique solution of (|6 
correct annular analyticity 
to ©. 

Let us next note that for the exact solutions of (|5l|)- 
(p3|), as Re£ — > +cxd, vanishes for the flat solution and 
has, for the other two 1/2 finger solutions 



v — > 1/A as Rc£ — » +oo, 



so that solutions to (63) with h regular as Re£ 

1&(A) - -p(l/A), 
for the 1/2 finger solutions, but 

= -p(Q) => * = 



(64) 
-co, 
(65) 

(66) 



For the analogs of the infinite fluid 

u n /n, * = AV>'(A) = e 2ni f, and 



for the flat solution, 
problem, with p = 
so nothing is changed for the flat solution's perturba- 
tions, but now for perturbations about the 1/2 finger, 
impedance (flux) grows exponentially in proportion to 
the instability's strength, and so are not free, and con- 
trollable by the pump that drives the flow: It is exactly 
as potent to control the change of flux as to control the 
growth of these instabilities. This is totally different from 
the infinite fluid case where these instabilities were freely 
impressible and uncontrollable. At this point we already 
mathematically conclude that the geometry of infinite 
fluid is physically incorrect, and the limit of large length 
clearly singular. This is the first main point we set out 
to establish. 

Let me emphasize this. The literature has been cal- 
culated in a physically incorrect geometry, and in conse- 
quence, all its conclusions lie suspect. (Its incorrectness 
in grander ways will transpire.) 

As for h's that are regular as Re£ — > — oo, all u's of the 
exact solutions diverge, and so the p's that determine 
them have constant limits (for simple powers, zero), and 
so these all have $ = 0. These will all turn out, as do 
the powers which relax as e~ 2nv , as relaxing stabilities. 
Thus, the flat solution is uncontrollably unstable, but 
the 1/2 finger solution can surely, under appropriate ex- 
perimental control, be reduced to non-exponential slow 
modes. It is our goal to show that these putative slow 
modes are non-existent. 

At this point we drop any further focus upon the flat 
solution, and restrict ourselves to the dynamical 1/2 fin- 
ger of (p3|), since ( |5^ ) is simply a long-time, long-fluid 
limit of this. We will now write down all possible "per- 
turbations" / upon this solution. 

Denote by x, y: 



= P c-c(A) , = „r_ 



x = y{-(, A) 



xy 



= e~ 2C = 1 + A (67) 
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so that 

v={l + x)/\ £=(l + y)/\. 
Denote by /3(A): 



/3(A) 



-2? 



1 + A 
1 + fcA 



by (3]38|). Finally, define 

i>(x,A) = /i(C, A), 
so that (B3|) reads 



and 



v(/3x, A) - v(x, A) = P(^y^) + ^(A), 



a(i/,£) = A) - A). 



(68) 



(69) 



(70) 



(71) 



(72) 



By straightforward differentiation, suppressing the A- 
dependence, 



1 + y 



f = (l + x)v'(x)- 
((1 + x)v'{x) + Xvx(x) - (x -> j/)). 



(73) 



The flux-bearing solutions, x — » 0, are, by summing (|7 
at the contractive (/3 < 1) fixed point at 0, 



Ud (x,A) = -5>(I±^)-p(i)); 
o 



A 



-A' 



(74) 



1> = -P(\) 



x pi x ] 



which we call "down-summing" , or "d£" . 

The flux-free stabilities, x — > oo, are, by summing the 
inverse, 

Vu ( X) \) = -J2(p(—^ )-p(oo)), * = (75) 



which we call "up-summing" , or "mS" . (It is immedi- 
ate to verify that the £ derivative of Vd + v u is precisely 
doubly-periodic when both sums exist.) Should p be nei- 
ther up nor down summable, then we simply split it into 
p — Pu + pd ~ V'p(A) with ip p , a piece of the ip for this 
p. (This is just a Laurent expansion when not explicitly 
possible.) Thus, we can determine all possible v's. 
Let us now write 



R(u) = -vp\v). 
For d£ solutions, by @, 

*(A) = -R(l/X). 



(76) 



To determine /, we need one last object, which by ( |69| ) 
is 



x J = a{x) = TTx- 



(78) 



Then, by (73), 



y ^ 1 - /3" + nx/3"a(A) 1 + /3"x 1 

TT^ TT^ a ~ R{ x )] 

This equation, in conjunction with ( |77| ) is remarkable: 
It says that the pump's measurement of extra requested 
flux as a function of time uniquely determines and is 
uniquely determined by the entire spatial shape of the 
flow: The pump can literally be programmed to create 
or reciprocally control any desired unstable conformation 
of the fluid! 

Next, in the finger regime (3 ~ \/kX and is astronom- 
ically, exponentially small, so that the sums in (^) are 
astronomically small, so that fd ~ -R(^), just the re- 
sult of the infinite fluid with one interface. For example, 
R = ln(l + av) produces a pole-dynamics singularity 



# = ln(l+-(l 
A 



with 



))=ln(l + -)+ln(l 



(80) 



1 + A/a 

which for a G (0, 2) is always £ s < for all A, and for 
a = 1 is just = C- Clearly, as (3 — > for a long 
enough channel during the finger regime, we have then a 
propagating finger of width (1 — e)/2, and all fingers are 
allowed. So it would seem. (It is false.) 

In any case, with R(v) analytic as Re£ — > +oo, these 
down sums are perturbations exponentially stronger (by 



roughly 



") on the driven interface than on the 



efflux interface, and clearly by ( |77| ) all instabilities. 

Finally, let us notice where all of /d's singularities are. 
First, R(Q) — since by (76) this is true unless p = 
lnzv + ... which by ( |55] ) means it contains a piece that 
changes the area of the fluid, which is illegal. Thus the 
poles at 1 -I- (3 n x — and 1 + j3 n y — are removable. The 
pole at y + 1 — 0, (or £ = 0) is a higher order correction 
to the basic singularity of / = ln£, and is of no special 
concern. Since (3 n — e~ 2 ™^, we then see that fd of (|7£ 
has singularities at 



(77) C+n = C. - 2n£ n = 0.. and C-„ = 2< - & n = 1.. (81) 



when -R(i^) has a singularity at Co = Cs- Most important 
is that C-o = — Cs is not a singularity. (By ([74]) the term 
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in the sum yields a = p{^)~P( v> ) an d / = vd v a dispatches 
the /?(£) term which contains the putative singularity at 
— Cs-) ln consequence — ( s can lie in physical fluid (0, £), 
and so, fd's can have singularities below, but arbitrarily 
close to ReC = 0, and so able to significantly modify the 
shape of the driven interface. Again, so it would seem. 

Also, in particular, fd is singular at 2£ — £ s if it is at 
£ s . This will prove to be of the utmost consequence. 

The up-sum solutions are of a totally different charac- 
ter. By (fra) we have 



fu 



P n + x 



A 



(82) 



—Y 



l-0 n + nxa(X) 1 + 0~ n x 
WTx ( A } 



It is useful to note that the sum of the two n = 1 terms 
in x produce 



U = -R{ 



l + (3- 
X 



(83) 



and f u in the finger regime is well approximated by just 
this leading term. As for singularities, with R{v) pos- 
sessing a singularity at ( s , f u is singular at 



C±„ = ±(C. + 2n0 n = 1, 2, .. 



(84) 



However, apart from £+1, by (|78|) all other terms are also 
singular at 



C±n = ±(-C + 2<) n = 1, 2, ... but not C+i- 



(85) 



Should R be logarithmic, (i.e. have cuts) such as (80) 
then the first term, (B3h is 



ln(l 



(86) 



This has two immediate consequences. First, £ s > — £ 
so that C s + 2£ > £, and so outside of physical fluid. 
(For £ s more negative than — 2£, one of the £+„ will still 
be in fluid.) Thus, such perturbations can only be im- 
posed after a time with A sufficiently small. Second, with 
Cs + 2£ > £, the only branching of (|86| ) that is reflection 
symmetric is with the cut running off to the left, and 
so throughout physical fluid, thereby violating channel 
width. Accordingly, branched i?(^)'s must contain paired 
singularities, with cuts beginning on one singularity, and 
terminating on a mating one (or ones). That is, if 



then 



R = Sa„ ln(l + a n v), 



Sq = 0. 



With ( p7|) enforced, and with R logarithms, the long fluid 
behavior of (|83f) is that of pole-dynamics with singulari- 
ties to the right of fluid (low pressure) sinks, but sinking 
zero net flux. 

For both second order perturbation and the theory 
of perturbations about perturbed solutions, we require 
some facts about first order solutions which follow by 
differentiating © and @. Differentiating @ by £<%, 
and setting v — > yields by (|l4|) and ( flij ) 

— / / (C)+/ / (-0=/'(2|-C) ! (88) 

7 

a result just as well obtained by differentiating ln([fl)|) on 

c 

Differentiating ( [45]) by £c%, together with a's antisym- 
metry produces, after noticing that 

= jtdt\(&(v), V ) = -/'(eiM^) 

and setting v — » z//7(£), 

—/(C) - /(-C) + /(2| - C) + 2*(A)/'(2£ - C) 

7 

= £//(£) = -#(£). (89) 
As a last observation important for second order theory, 

by ©, 

P = *(A)/' = -i/fl^ and P* ee *(A)/ W = -£A^ ( 90 ) 
so that 

£<%P = vduP*. (91) 

Notice also that 

vdpu\ c _t = vd u u(-^—,t;) 



7(0' 



and 



= (£^ + 

6 4 t 7 7 (£) 

It is now straightforward to write out second order 
perturbation theory. First calculate the pull-back of G 

of © to c - i By m 



l(C-£)=9(-Q = U + P)(y,Q, or 
7 



<7M = (/ + P)K(£),£). 

For <7* = g(£, v), we have 

s'(C-D - </(£, -) = (/+P)K7(r). -) = (/+-P)(2£-0- 



(87) 



7 



7 7 
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Now, using ( |89| ) and (Bsj) multiplied by obtain 

iU'(C-l) = 



(/ + P)(f - P l ) - ¥-(f + p f - Uf + p ) R (0- 



27 



Combined with the first term of ( pq ) pulled back to £ — £, 
we have 

g(c - - (/ + + ^p*)+ 



»U\U + P)(f -p*- R(0) + 

We now equate this expression to G(( — as obtained 
from f p , (|29]), using ( |l5| ) multiplied by ^. Defining the 
second order flux 



(92) 



after obvious cancellations, we produce an expression of 
the form vd v { ), and so, with an integration constant 
purely of £, obtain the fully integrated second order 
analog to ( |38| ) 

A(v,£)-ip{\) = (93) 

i/(P* - R(0) + \p{f - R(0) + A g (l,0). 
z z 7 

We now proceed exactly as we did to achieve ([30]), with 
now an extra right-hand driving term, ip and fj, are chosen 
as the minimal required terms to produce an A analytic 
over physical fluid. Here we must generally avail our- 
selves of an appropriate doubly periodic homogeneous 
solution, should the particular solution determined by 
the driving term from first order happen to produce un- 
acceptable singularities. (The form of ( p3[ ) with /*— i£(£), 
and the construction of F from the potential A + 5//* 
assists in eliminating such problems.) 

To end this section on formalities, let us sec how 7, or 
r is modified to first order in e from 7 and T. Expanding 

/(-C)-/(2£-C)-2e*/'(2£-0 
+e(/(-C)-/(2|-0) 
= r(/(C) + £ /(C)) + 0(e 2 ). (94) 

With /(-C) - /(2| - C) = f (/(C)), and dividing by e, 



/(-C) - /(2£ - C) = -(T(/(0 + e/(0) - r(/(0)) 

e 



+2*/'(2e-C) + 0(e). 



r'(/(0)=^(C), (95) 

7 

by (p9|), the perturbation theory, 

- e (nAO+ef(0)-Hf(0)) = R(ef)+t'(f(0)f(0+O(e). 
Taking the limit e ~ (r must be analytic), 

T(f)=t(f) + eR(e f ) + 0(e 2 ). 

That is, for / p to be the full perturbed solution to our 
problem, and with £ = e^ p ^\ 

r(/p(C)) = f (/ p (0) + eZ2(0 + 0{e 2 ) (96) 



7(0= 7(0(1 + eR(£)) + 0(e 2 ). 



(97) 



This result has deep consequences for perturbations, 
which we now spell out. 



5. DERIVATIVES VERSUS PERTURBATIONS 

The idea of perturbation theory is that given a so- 
lution, wc find all the (partial in e) derivatives of the 
equations of motion about this point which, in principle, 
determines a power series for the solution with leading 
term the known solution. Should the solution actually 
be differentiable over some compact region, then the first 
order solution / + ef, (/ just a particular tangent) with 
e small but finite, agrees with an actual solution over 
the region within a uniform 0(e 2 ) error. In general one 
needs either powerful bounds or the entire series to all 
orders with the knowledge of its convergence, to know if 
this derivative actually exists. It is only then that / + ef 
is an approximation to a solution. And this can almost 
never be settled to any finite order in perturbation with- 
out exceptional extra knowledge. The exact result (4j3^) 
and it uniquely determined perturbed form, (4.|96|) pro- 
vides us with some such knowledge. 

Specifically, (4j97j) says that if / + ef for a finite, fixed 
value of e is an approximation to an actual solution, then 
7 exists as a function purely of £, and the leading term 
and part of its O(e) corrections resides in 7(0- That is, 
in (3.193) what must be true is that the first term is 



and not 



r(/ + ef) 



f (/) + ef '(/)/ + 0(e 2 ). 



(1) 



(2) 



Since, by (35), 



Evidently, if f is always bounded over the region we care 
about, then (||) is altogether equivalent to (|l|). Should / 
not be bounded over the region of control then only (jl]) 
applies, and (0) is wrong. The formal 0(e) perturbative 
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result, (4.|9|), relies upon (§), not ([[]). That is, if the per- 
turbation / entails singularities there are deep potential 
problems unless e is viewed as an actual infinitesimal. 

Now, what we might have liked to say is that the re- 
gion of concern is just [0, £], the domain of physical fluid, 
so that even if / has singularities outside this region, 
we need pay attention only to a domain free of them in 
which (||) and (|J) are equivalent. That both interfaces 
must be satisfied, with its specific form of the method 
of images, (4j3^) prevents us from determining satisfac- 
tion just within [0, £]. By the equations of motion for 
the driven interface, singularities within [— £, 0] must be 
compatible (e.g. have residues correctly determined to 
0(e)) with the analytic form of f p within [0, £]. But also, 
the simultaneous satisfaction of the g equation, yield- 
ing (4.pS6|) requires including the behavior of / up to a 
least 2£, if not even 3£ (since — £, near a singularity, can 
approach £). Thus, f p ~ / + ef must consistently agree 
with singularity structure over a region basically [— £, 3£]. 
Indeed there is every reason to expect f + ef to arbitrar- 
ily disagree outside this domain, so far beyond physical 
fluid. 

Let us first see that the up-summed solution, stabil- 
ities with their far-away singularities, is in conformity 
with (4.^) and (]l|) throughout this domain, even at sin- 
gularities. Within this domain, only the leading term 
(4.[S3"1) possesses a singularity, at £ s + 2£ when R(y) has 
a singularity at ( s . Since in this case Q s S [— £, 0] in the 
worst case (strongest perturbation), f p must be analytic 
at — ( s , which is then in physical fluid. So, in (4.|3^) we 
consider 

C = -£ + X, 
f p (Cs -x)- f P (2t; g +( s -x)~ T(f p (-Cs + x)) 

~t(f p (-t s + x)) + eR(eK-S>+^). (3) 

Now, f p s perturbing part is singular at 2£ ff + £ s by 

(4.|||), f p is regular at ±( s , as is Q7p(— C*)))j but R 
is also singular at — ( s . Indeed, (4j83j) of perturbation 
theory exactly agrees with this, and so, we conclude that 
up-summing, even with singularities, is compatible with 
f p existing as an approximate solution. 

Next we turn to down-summing. With R(v) = v n , n > 
0, it is elementary to exactly sum (4j7|), while by (4^), 
such a perturbation is flux-bearing with $„ 



-A" 

rapidly exponentially growing, but visible to the 
pump driving the flow, and so in principle, controllable. 
These perturbations are nonsingular, and so for small 
enough e (|j) and (||) are equivalent. At the end-points of 
the domain we must consider, the perturbing /'s become 
enormous, so that f + ef is probably reliable only for 
e ~ e~ L . For larger e, when the perturbation has grown 
significantly with ip, we have no idea what its nonlinear 
fate is. 

Finally, we consider summing the i?„'s into singular- 
ities, to produce, for example, pole-dynamical R's such 



as ln(l + av). Let us now see that if these singularities 
lie within [— £,0], then the perturbative result is not an 
approximation to any solution of the problem. The test 
is for 

C = Cs + X 

where, for C, s G [— £, 0], — Cs is within physical fluid. Here, 
(4.||) is 

f p (-( s -x)- f p (2£ g -Cs- x) = T(f p (Cs + x)), 
or, with (4j9q), and picking off singular terms, 



- e/ d (2£ - C s - x) ~ T{f{Q + ef d (( s + x)) + eR(f(Q). 

(4) 

Now, R is regular at /(Cs) (it's singular at /(—£«))) but 
F is to be evaluated at an infinite value of its argument, 
at which, in fact, it is regular. However, by (4.|7^), 2£ — ( s 
is a singularity of fd, and indeed, within the domain of 
concern, only Cs and 2£- — £ s are singularities. But then, 
(^) is inconsistent. (Had we used the incorrect form, (Q), 
the two sides would have agreed with identical residues.) 
That is, down- summed solutions are compatible with be- 
ing approximations provided that R(v) is free of singu- 
larities in [—£,0]. This implies isolated pattern selec- 
tion. To produce a finger of different width, the leading 
pole-dynamic term R{y) = ln(l + av) is required for val- 
ues of a not exponentially small with fluid length L. In 
this case during the finger regime A < 1, M > 1, or 
1 <C tp <C L, the singularity ( s must move 0(e~ 2v ) near 
to ( = for the logarithmic stretching to take place that 
modifies finger width. To say that f + eR for fixed finite 
e is an approximate solution is that 

f p ~ / + eln(l + a/A) + eln(l - e^), 

so that near £ s , 

-< 



f 



so that 



fp,v f S _ fi 

r, S SS Ss' 

Dividing the equation of motion for the driven interface 
by f p (C), we have, just as for pole dynamics 

f P (-( s (<p),v) =x s +0(e 2 ), x' s =0. 

That is, if ln(l + av) is to produce a changed finger, f p 
near Re£ = must include the singularity and f p (— Cs) 
must produce a stagnation point. f p near its other near- 
est singularity, presumably at 2£ — ^ s must be compatible 
with the method of images (4.p6|). That is, (4.|3q) must 
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be satisfied for |Re£| < S with S — > + . This is exactly 
what we have demonstrated is false in (Q). 

Indeed, with singularities never within [— £, 0], the per- 
turbation is exponentially small with fluid length during 
the finger regime, and so for a long enough body of fluid, 
without impact on the solution /. That is, we have ruled 
out the possibility raised at the beginning of Section 4. 
that requesting an exactly translating solution might be 
too brittle a question to determine selection for a finite 
body of fluid. Furthermore, with / warding off singu- 
larities in [—£,0], or roughly by the length of physical 
fluid, the greatest host of perturbations that could in- 
duce finite-time singularities have also been warded off. 
All that is significantly left are the exponentially grow- 
ing v n perturbations. Left to themselves, these rapidly 
lead / into a nonlinear regime about which we know lit- 
tle, but surmise tip-splitting. On the other hand, as they 
proportionately modify flux, they are controllable. It is 
conceivable that a simulation of this flow constrained to 
prevent rapid flux variations might, near enough to /, be 
stable. 

Apart from these exponentially growing fluxes, there 
are no slow modes to modify £ ~ Lq — ip to £ ~ 
Lq — {1/\— 1)^. That is, this flow protects the far down- 
stream behavior of the 1/2 finger without surface tension 
on the efflux. The ability of / to ward off nearby singu- 
larities lies in its ability to satisfy (4.|36|) without any im- 
age singularities present, whereas all putative nearby per- 
turbations with singularities must possess them, and are 
then inconsistent. Clearly, with surface tension present, 
a method of images reminiscent of (4.|36|) must exist, and 
so one can surmise that the appropriate S-T solution is 
again special. Even though surface tension is a perturba- 
tion that spoils the critical symmetry of this paper, there 
is now a reinforced reason to surmise that (3.pT|), 

1 - 2A = (2n) 2 B e , 

at least for A near 1/2 is correct, which I offer as a con- 
jecture. 

6. PRECIS AND DISCUSSION 

We have followed a clue laid down by reflection sym- 
metry that the analytically continued equations of mo- 
tion are symmetric under parity, £ — > — Q. This suggests 
that far upstream (Re£ < 0) singularities, which deter- 
mine the shape of the driven interface at Re£ = 0, are 
related to and co-determined with far downstream be- 
havior, which is the nature of the efflux in a physical 
experiment. To see if any such thing holds, it was then 
mandatory to consider the flow of a finite body of fluid. 
The simplest feasible such problem, fully posed within 
the purview of the analytically continued conformal ma- 
chinery, consists in replacing the Riemann mapping of a 
disc to an annulus. The case of the simpler disc embraces 



all prior studies, and corresponds to fluid going off in- 
finitely far to the right with a simple pole at +oo. In the 
annular version, the available channel is still infinite, but 
fluid terminates at a downstream interface, the behavior 
of which constitutes the efflux for this configuration. 

Requesting a purely translating driven interface re- 
ceives a unique solution in the finite problem: There is 
sharp, precise selection for this question. Examined in 
pole-dynamics, even including singularities to the right 
of the efflux interface, precisely and only this selected so- 
lution exists. Examining a class of solutions that exist 
for the problem of the disc, but not here; namely the 
usual Saffman- Taylor solutions, we notice that selection 
is expressed on the efflux interface, even when asymptot- 
ically far downstream, with only the 1/2 solution main- 
taining it at a fixed pressure. Solutions with A < 1/2 
also correspond to solutions to the finite problem, but 
with surface tension now smoothing the efflux interface 
(although, of course, with zero surface tension still on the 
driven interface). The explicit relationship, 1 — 2A = a jv 
provides the physical interpretation of all the S-T solu- 
tions. Moreover, that the infinite fluid version of the disc 
allows all these possibilities with equanimity, determines 
that its pole at +00 is truly within physical fluid, so that 
the problem of the disc is the theory of an unfinished, 
unterminated experiment, and is not the limit of a very 
long body of fluid. That is, the physical limit L — > 00 is 
singular. 

We end up showing that solutions to the finite problem 
near the selected one all preserve the efflux condition 
1 — 2X = = a/v. Although we have not worked out 
the theory for surface tension on the efflux, considering 
that the A / 1/2 solutions do meet this condition, and 
that the condition is conserved under perturbations in 
the A = 1/2 case, it is the most natural surmise that the 
width-efflux condition is generally correct, at least for A 
not too much below 1/2. 

Most importantly, we reach the conclusion that while 
many consequences of the distant efflux are exponentially 
small, this is false for pattern selection. This is so differ- 
ent from the generally presumed behavior that for this 
reason alone, the annular configuration we discuss should 
be put to the experimental test. It is important in this 
regard to note (private communications from several of 
the experimenters) that this selection with beautifully 
formed fingers is not so easy to experimentally produce, 
generally requiring fiddling with the efflux to become re- 
liably reproducible. For example, Tabeling et al || suc- 
ceeded only after fitting an "impedance-matching" plug 
to the efflux. It is regrettable that the experimenters 
haven't carefully observed and calibrated such dependen- 
cies. This is not for any failure of the experimenters, but 
rather, illuminates an important aspect of the relation- 
ship between theory and experiment. Namely, unless one 
has prior insights - basically theoretical knowledge - then 
one doesn't quite know what is worthy to observe. Nec- 
essarily, in order to better observe what one has chosen 
to, the experimenter exercises his art to control what ob- 
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fuscates. But some of these thrown-out details might be 
at the actual heart of the system being observed. How 
is one to know in advance? And certainly, until certi- 
fied by experiment, there is no way to be sure that the 
theory considered is significantly related to the physical 
phenomenon. 

Having determined that at least for the elementary S-T 
solutions that only A = 1/2 has any connection to a phys- 
ically finite flow with efflux at atmospheric pressure and 
no surface tension on it, we go on to question if asking for 
pure translation in a finite problem is perhaps too brittle 
a request. In the flow of a finite piece of fluid, there is a 
finite, definite period of time over which the solution is 
that of a finger with a large body of fluid still downstream 
from it. The driven interface at velocity 1/A penetrates 
into the fluid whose terminating interface, virtually flat, 
moves with velocity 1. For A = 1/2, by the time the 
fluid has translated its resting length L under the mean 
velocity of its flux (i.e. 1), there remains but a thin layer 
of fluid still downstream from the driven interface, as 
the efflux interface begins to curve into a similar finger. 
Calling this finite period of time when the pure finger 
is well-formed (i.e. with long horizontal sides) but still 
with a large body of fluid behind it the "finger regime" 
(as opposed to afterwards, the "breakthrough regime"), 
we can ask if there are solutions that within measur- 
able error appear as perfectly translating, but actually 
possess deviations exponentially small with L, only be- 
coming physically important in the breakthrough regime. 
Should this be the case, then the strong selection is a mis- 
leading artifact of too mathematically posed a question. 

Since there are no exact solutions we know of (e.g. 
all versions of pole-dynamics) save for the pure 1/2 S- 
T finger, we regrettably had no choice but to perform 
a perturbativc analysis. As we erect this machinery, we 
are immediately led to notice that exceptionally useful 
new variables, basically (/(C, <p), /(— C> <p))> constructed 
from the old ((,ip), are guaranteed of invertiblity, since 
the equation of motion for the driven interface is pre- 
cisely the statement that the Jacobian is the constant 2. 
Translating in ( an identical fact for the equation of mo- 
tion for the efflux interface, we notice that given the first 
equation of motion, the second can be fully integrated 
into the relationship 

/(-C, ip) - m g - c, <p) = r(/(c, ip)). (r) 

What is crucial here is that V depends only upon the 
one new variable /(£, ip). This formula serves as the pre- 
cise way in which the method of images is to provide a 
solution compatible with both interfaces. 

In these new variables (more conveniently, their expo- 
nentials, (£, v)) perturbation theory is especially simple 
for just one interface. Together with (r) expressed per- 
turbatively, we are able to fully integrate the perturbed 
p.d.e.'s for both interfaces, and reduce the problem to a 
purely algebraic one. We explicitly did so here for the 
first two orders of perturbation, and to first order wrote 



down all possible solutions. The most striking such solu- 
tions are just the usual infinite fluid pole dynamical ones, 
but decorated with exponentially small image terms, van- 
ishing with L — > oo. Since we already rigorously know 
that there are no such solutions, we wonder how they 
perturbatively appear so naturally and fluently. Before 
discussing this conundrum, we notice something central 
to the heart of the finite problem. 

A finite body of viscous fluid has the natural con- 
stituentive relationship of impedance, the ratio of the 
net pressure across the fluid to the flux of its transport. 
Although impedance is infinite in the usual infinite fluid 
problem, it is less clear that its time derivative should also 
exactly vanish, as it does in the theoretical literature, but 
ostensibly in disagreement with experiment. This is a sig- 
nature that the infinite boundary geometry is unphysical, 
even construed as the limit L — > oo. Should a deforma- 
tion of the fluid's flow produce a change in impedance, 
that is to say a change in the conserved flux through the 
entire fluid when pressure across it is maintained fixed, 
then such a deformation requests power from the pump 
that maintains this pressure difference. That is, such 
a deformation is driven by the pump, is controllable by 
the pump, and is certainly not spontaneously impressible. 
Such a circumstance is hardly surprising since a deforma- 
tion is not just on the driven interface, but throughout 
the incompressible fluid that must everywhere adjust to 
accommodate to a distortion of its boundary. This is 
very different in the infinite fluid geometry, where the 
promiscuous pole at infinity simply picks up whatever 
need be, since for the Riemann disc, the disposition of 
the one interface uniquely determines everything. 

Performing the perturbative analysis of the finite fluid, 
we are forced to immediately decide if the perturbation 
is to modify flux, and so be externally driven. Should 
this not be allowed, we discover that the unperturbed 
finger has only purely stable modes, although not enough 
to allow for the arbitrary independent perturbations of 
both interfaces. The residual modes, all instabilities, are 
all driven by the external pump, and so, in no sense free 
fluctuations. These same modes, as L — > oo, are precisely 
all the modes of the infinite fluid problem, although in 
that configuration they are all free. This now cements the 
observation that L — > oo is singular, and the infinite fluid 
geometry a physically wrong limit. There is no reason to 
doubt that these exponentially growing modes are real - 
they are, however, controllable. Indeed, a record of the 
flux through, or the power delivered by the pump over 
the course of time uniquely determines the disposition of 
this unstable flow throughout the entire body of fluid, 
and throughout this interval of time. 

Granting these modes and an experimental control 
that allows them to grow at most algebraically in time, 
we now ask about the slow modes produced by summing 
all these fast modes. We especially care about summing 
them into logarithms, which then modifies the linear 
change of impedance in time, which then modifies the ef- 
flux transport from L — tp to L—(l/X—l)<p, hence modifies 
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A from 1/2, and so, with no surface tension on the efflux 
produces a flow that violates 1 — 2A = a /v. Here, resort- 
ing to (r), we discover that these logarithmic modes are 
in fact purely formal, where with e the strength of pertur- 
bation, they fail to be solutions for any finite value of e, 
no matter how small, thus resolving the conundrum that 
they don't exist in any correct non-perturbative analy- 
sis. With this, we establish that the solution of A = 1/2 
is definite and isolated in the neighborhood of the pure 
unperturbed finger, and that we had not posed too brit- 
tle a question in asking for exact translation. Evidently, 
one would far prefer a purely non-perturbative treatment 
for a problem with a rather subtle resolution. Neverthe- 
less, this work serves as a tonic against the use of linear 
response methods in such delicate systems. 



7. APPENDIX: SELECTION WITHIN POLE 
DYNAMICS 



The ^-translation invariant solutions of (3.37), which 
we shall here refer to as /„ with n as in that equation, be- 
long to a certain well-studied class of solutions (Class (27) 
of Reference Q ) We used no properties of that Class, 
but rather deduced (3]37]) from dynamical symmetry ar- 
guments. This Class has a most important property: its 
members are exact solutions of (2.14'). These particu- 
lar dynamical solutions are the so-called "pole dynam- 
ics" . We shall, to render this paper more self-contained, 
re-derive these results here, as they are so trivial un- 
der reflection-symmetry, but then easily extend them to 
more complex variants. Most importantly, we shall see 
that even within this significantly enlarged class of pole 
dynamics, the only solutions for both interfaces within 
the class are precisely the f's of (3.[37]). That is, there 
exist no perturbations to / whatsoever in this general 
class, which is to say that had the dynamics been ex- 
hausted by pole dynamics, then / is rigidly isolated as 
the solution, and hence sharp pattern selection. 

We start with the delineation of pole dynamics 
(Class(27)): 



/ = /%) + C + Sa k ln(l-e Cfc - c ); 

ReCfe(y) < 0; oifc = const 



(1) 



where / must satisfy (2.14') and be reflection-symmetric 
(P real; (ak, Cfe), (<5fe, Cfe) both present). Paying attention 
to (2.14'), write 



/' = ! + £- 



.C-C* - 1 



, U = & - s - 



>C-G 



(2) 



Ca 



Consider C - Cfe- Then 1//'(C) - 0, f v /f'(Q 
-C k (<p). Dividing (2.14') by /'(O/'K) and taking ( 







-C k + U(-(k)/f'(-(k); 



o = U(-Ck,<p)-Ck(<p)f'(-Ck,<p) = -rzfi-M 



dip' 



f(-(k,v) = z k= const. 



(3) 



These points, z& = /(— Cki'fi) are "stagnation" points, 
in the sense that as <p — > oo all Re Cfe — > _ , so that 
Zk is very near to a point on the interface C = is , and 
the interface may never pass through this point. More 
correctly, the interface must asymptotically come to rest 
at 



z k - a k In 2 



(4) 



since as Re Cfe — > — Cfe — > Cfe > and / is singular at Cfe- 
Setting C = ii]k + it£k with Cfe = Cfe + iWki writing down 
/(C), /(-Cfe), and subtracting yields 



f[ir]k + it£k) ~ Zk - ak In 2 + ak ln(l - it). 



(5) 



For a flow within an initially almost flat interface that 
grows wrinkled, all the Zk have real parts to the right of 
the interface, and so, by (|J) must have Re ak > 0, since 
the interface can never have passed through Zk itself. 

(|J) determines all the Cki'fi) in terms of j3. f has deriva- 
tives analytic at Re C —* ±oo: By (||), 

/'(+«>) = 1, /'(-oo) = 1 - So, 

f ip (+oo) = (3\ / v ,(-oo)=/3' + Sa fe C^ (6) 

and so by (2.14'), 

(1 - Ea)j3 + (/? + Sa^Cfe) = 2<P + const, (7) 

thus fully determining /. 

As ip — > oo C(, — > 0, as Re Cfe — > 0~, and by f's imaging, 
Im Cfe — » 0, 7r only (see Paper I). By (||) outside arbitrarily 
small disks about the (k,f v ~ (3', and by (2.14'), 



/3' 



no+n-o 



with /' - 0, so /'(C, v) - /'(C). Hence 



1 



/?' = — = const, 
A 



(8) 



and /'(C) + f'(-C) - 2A, /(C) - /(-C) ~ 2AC and so on 

C = is, 



y(s) = Im/(is) ~ As s G (e, ir — e) 



(9) 



which is to say it is a finger of width A of the channel. 
By (§ and (|) 



A-i = i(l-£a) = i/'(-oo). 



(10) 



Now, the /„ of (3.^) are explicitly of form (|l]) with 
/'(— oo) = 0, and so A = 1/2. We ultimately care about 
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/i , a single 1/2 finger in a 2tt channel. It is now natural to 
ask if there are other Class (27) solutions for both inter- 
faces, since then we could analytically solve the problem 
of such perturbations to /. It is easy to see that there 
are none. 

First, substitute (2.[l|) into (2.14") to obtain 

2 - /'(-c + e 9 )/ v (c + e 9 ) - /'(c + o/ P (-c + 

= ^ g f'(t + Qf'(-t + Q- (ii) 

with £ g finite, taking the limit Re C ~* +oo, where 
Class(27) is analytic, the left hand side vanishes, as it 
is the limit of simultaneous (2.14'). Hence the only /'s 
that are allowed satisfy 



£»/'(+oo)/'(-oo) - 



(12) 



But /'(+oo) = 1, and for any non-flat solution ^ 0, 
and so /'(— oo) = 0, A = 1/2, and so, 



Ea = 1 



(13) 



is a consistency condition for / to obey both equations 
of motion. Thus, A = 1/2 is selected. 

There is an important simple fact related to why 
=/= 0. There is a unique conformal (i.e. invertible) 
map h from the fluid at time ip to a rectangle with sides 
£ 9 and 2n. The aspect ratio, £ 9 /27r, of this rectangle 
is termed the "module" of the region mapped, and is 
uniquely determined by that region. Moreover, it enjoys 
an exact estimate: 

2< s ( ¥ >) < A(<p) = A 

with A the area of the fluid, the constant in time Aq. 
Unless h! = 1, the inequality is strict. Thus, with A Q = 
2ttLq, we have 

That is, for other than the flat flow £ g (<p) is stricly smaller 
than the resting length of the fluid. In particular, for a 
solution which was flat in the far past, £' g ^ 0, since it 
had the value Lq in the far past, but certainly below it 
at finite times. 

With / of (H), g of (2.p|) is 



(/3 + G)+C + Sa fe ln(l-e 



(14) 



Between / and g, the £ — > oo equation, (Q) is identi- 
cal with Ea = 1. So, there is one such equation for f3. 
However together with (0), we now also have 



(15) 



With a real (k , by f's reflection symmetry such a (||) , ( |l5| ) 
is a real equation. For each complex each is complex, 
and hence a pair of real equations. So, with m complex 
£fc's and n real ones, we have 



# real eq's = 4m + 2n + 1 



(16) 



and 



#realvbl's = 2m + n + 2 (0/s + £ 9 + /3) (17) 
and so the system is overdetermined by 

excess =2m + n— 1. (18) 

So, save for the case of just one = C which is real, and 
even with just one complex £, the system is overdeter- 
mined. The one real case is just /i of (3]3S|). 

Let us now show that the only solutions that become 
flat as ip — > — oo are just the /„, with the overdetermina- 
tion in fact inconsistent in all other cases. 

Writing down (||) for a and the difference of jl5| ) 
and (^), we have 

z k =/9-a + Sa^ln(l-e Ci+Cfc ) (19) 

/I _ e -2£ s +C*+C< 
z fc = A fc = 2£„ + Hat In ' 



For all Re (k — > —oojj. 
immediately have by (EC 



A k = 2£ . 



o°),S,g -> £o (finite), we 
(21) 



Expanding (gOj) , since all Re < and £ g > and so 
the series are absolutely convergent, produces 

oo 

Ea £ ^i(l- e - 2 <")e'' (a+Cf) 
which can be resummed as 



E 



-(l-e ar6 «)(Eo<e rC£ ) 



So, 



2(^o-e 9 ) = E u '-^) era ; 



-ft?) = -(l-e- 2r «»)Ea^ 



(22) 



The first of fl22j) implies that 

u r [ip) = r < n; u n ^ 
and e™^* = e n< >° for some n > 1. 



(23) 



Thus 



Cfc = Co + 27ri- k = 0, ....,n — 1. 



(24) 



For u r <„ = 0, by the second of (p2|), J2"=o a i^ % ™ 
, r = 1, n — 1, and so by finite Fourier, with (|l3|), 

a fc = 1/n. 
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Finally, taking Re C sufficiently large and analytically 
continuing the result, 



1 n— 1 

/ = a + c + - y in(i 



e Co-C e 27rifc/™) 



fc=0 



oc / 1 n— 1 

/3 + C - V ie r «°-« - Y e 2mkr ' n 
^— ' r \ n ' 

r=l \ fc=0 



1 1 



n ^ £ 

=i 



/3 + C+-ln(l-e n(c °- c) ). 



(25) 



By reflection-symmetry e™^° is real, or e Tl ^° = ±e ra ^ for C 
real, and so 

/ = /3 + C + - ln(e n(c -^ ± 1) 
n 

which then to satisfy (2.14') is just /„ of (3.|37|). It must 
be remembered that (2.14') is nonlinear: the only /'s of 
Class (27) is precisely one of these /„'s not any linear 
combinations whatsoever. That is, there can be no Class 
(27) perturbations to /, so that within this Class, / is 
perfectly stable of unique 1/2 width. 

Now, Class (27) arose in discussions of the infinite 
channel: in the natural variable e - ^ = ui, u = is 
analytic at uo — 0, where it has a simple zero. e~f is the 
unique (under reflection symmetry) Riemann map tak- 
ing u = to lo = of the physical region to the unit 
disk. With two free boundaries, the map is to the annu- 
lus e~^ s < \u>\ < 1, and so is generally not analytic in 
\u\ < e-ts. The fs of (§ are a Taylor expansion about 
uj = 0. For two free boundaries, it is generally to be a 
Taylor-Laurent expansion. Let us write this as 

f = /3 + C + Sa fe ln(l-e Cfc - c )-Sa+ln(l-e-^ +c ) 
with ReCfe < 0, ReC^ > 0. (26) 

To fully consider these solutions, we extend the pole- 
dynamics of (^) to include the special circumstance that 
the Cfc and might be matched with = —Cfc, which 
is the "closest" that ( p^ ) can come to a solution (There 
are none.) 

Writing the derivatives, 



/' = 1 + E 



e<-& -1 ' e tf-C-L 



U - - s- 



a k (k „ a k C k 



(27) 



First 



/'(+oo) = 1 - Ea+, f(-oo) = 1 - Ea, 
U(+co) ={3' + Ea+C,+', U(-oo) = P' + Za k (' k , (28) 



so that (2.14') yields 

(1 - Ea)G8 + Ea+C fc + ) + (1 - Ea+)(/3 + Za k ( k ) 



2ip- 



const. 



(29) 



It is next immediately clear that (^|) is unmodified for 
both Cfc and */they are unpaired: 

/(-Cfc, <p) = z k , f(C+, <p) = z+, C fc + C fe + + 0. (30) 

Should Cfc + Cfe~ — (arrange the indices so that mat- 
ing C's have the same index), then the derivation of (|^) 
fails. Instead we compute with ( p7j ) /'(±(Cfc + e)) and 
/^(iCCfe + e ))> and proceeding with a little care deter- 
mine that 



Zk = a k f R (-( k ,(p) + a£f R (( k ,(p) 



(31) 



where f R means drop the term that diverges at the spec- 
ified argument of / precisely of form (|26|). Most impor- 
tantly, there is just one equation ( pl[ ) for the matched 
pair ±Cfc which lowers the overdetermination count of 

The problem is that it doesn't lower the count enough, 
since if a pair matches in fit can't in g and visa versa: 
The singularities of g, by (g6|) are at Cfc — £, g and C^ — £g, 
so that if Cfc and Cfc" are paired, and hence with mean 0, 
in g they have mean — £ ff and are no longer matched, and 
so each produces its own equation (|i"5|). This, together 
with the notion of the "method of images" , determines 
the best possible way to place /'s singularities. Consider 
a real 



Co = C < o. 



(32) 



With — C > and approaching as <p — > +00, it is almost 
impossible for — C not to be in physical fluid, and so take 
it as unmatched. Next consider more Cfc: 



C fc = C - 2fcC fl 



0, 



.,71. 



But then 



-Cfc 



-C + 2fc£ 9 fc=l, 



(33) 



(34) 



are all matched, greater than 2£ ff , so certainly outside 
physical fluid. These are fs singularities. For g, we have 



Cfc 9 = Cfc-i-e s = C + 4-2fc4 k=l, 



,n+l (35) 



and 



/-+g 

Sfc 



Cfct-e, 



-C - C s + 2k£ g 



< 9 k k=l, 



(36) 



Again all the Cf are matched, save for Cn+i- Thus both 
/ and g have n matched pairs and one unmatched singu- 
larity. (This is the best that can be arranged within the 
£ ff annular domain.) Counting, this is 2n + 3 real equa- 
tions for 2n + 1 singularities, £ g and (3. This has a right 
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sense to it, save for the difficulty that the 2n+ 1 singular- 
ities are just combinations of £ and £ s , so that there are 
precisely always 3 variables, and so, overdetermined save 
for n = 0, which is simply f±. It should be reasonably 
clear that / of ( p6| ) affords no exact solutions other than 
the /„ of (^-translation invariance. That is, there is no 
pole dynamics at all. Rather, that class, natural to the 
dynamics, truly arises as a guessed- at extension of the 
Saffman- Taylor stationary solution, but exists only for 
1/2 fingers in consequence of displacement invariance. 

That is, the 1/2 finger has unlimited stability and se- 
lection within the scope of any conceivable pole dynam- 
ics. Insofar as the finite-time singularities in the litera- 
ture are all of pole dynamics, /i is evidently totally free 
from them. 
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